Ferromagnetism and Antiferromagnetism of Correlated Topological Insulator with 

Flat Band 
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In this paper, based on the mean field approach and random-phase-approximation, we study the 
magnetic properties of the spinful Haldane model on honeycomb lattice of topological flat band 
with on-site repulsive Coulomb interaction. We flnd that the antiferromagnetic (AF) order is more 
stable than ferromagnetic (FM) order at (or near) half-fllhng; while away from half-filling the phase 
diagram becomes complex: At large doping, FM order is more stable than AF order due to the 
flatness of band structure. In particular, we flnd that at quarter filling case, the system becomes a 
(5 = 1 topological insulator which is induced by the FM order. 



I. INTRODUCTION 

From 1980 year, topological ordered states become 
a hot topic and a lot of efforts had been done on 
this issue Integer and fractional quantum Hall ef- 
fect (IQHE and FQHE) have attracted great attention 
and been studied a great deal since first discovered in 
two-dimensional electron gas with Laudau levels induced 
by strong magnetic field]^, @|- Haldane's honeycomb 
modelQ demonstrates IQHE in lattice model without 
Landau levels while FQHE had not been discovered in 
lattice models until recently [^-Q- The topologically non- 
trivial flat-land was pointed to be key in realizing FQHE 
in lattice models based on the mechanism of quadratic 
band touching. And the lattice models are designed 
reaching a high flatness ratio (the ratio of the band 
gap over bandwidth) about 50 in different lattice models 
kagome lattice, square lattice and honeycomb lattice) [B- 
In these models, the hoppings of next-nearest neigh- 
bor (NNN) or next-next-nearest neighbor (NNNN) be- 
sides nearest neighbor (NN) are introduced to achieve 
the flat-bands. And the flat-bands of these models are 
shown to be with non-trivial topology (with non-zero 
Chcrn number). To obtain FQHE states, people always 
consider these tight-binding models of spinless fermions 
or bosons with the NN and the NNN repulsions under a 
hard-core condition. Then the numerical results confirm 
the existence of topologically nontrivial flat-band [S!]. 

On the other hand, the spin orders are important mag- 
netic properties in a strongly correlated electron sys - 
tem, on which people had also paid much attention jlOl|. 
For example, an interesting problem is about the ferro- 
magnetic (FM) order in correlated electron system with 
flat band. According to the Stoner criteria, for corre- 
lated electronic system with flat or almost flat band, 
the Coulomb interaction will remove the band degener- 
acy and drive the system into a ferromagnetic ground 
state [ni [Hi- So we may ask a question : in a spin- 
ful fermionic lattice model of topologically nontrivial 
flat-band with on-site Coulomb interaction, whether the 



ground state is FM order? This is the starting point of 
this paper. 

After doing mean field calculations, in a typical two- 
component fermionic lattice model of topologically non- 
trivial flat-band with on-site Coulomb interaction (the 
spinful Haldane model on honeycomb lattice [3. ITsl Il6j). 
we flnd that the AF order is more stable than FM order 
at (or near) half-filling; while away from half-filling the 
phase diagram becomes complex: At large doping, FM 
order is more stable than AF order due to the fiat band. 
In particular, we found that at quarter filling case, the 
system becomes a Q = 1 topolo gica l insulator which is 
induced by the FM order. In Ref.[l3|, it was proved that 
the ground state of the Hubbard model with topologi- 
cal fiat-bands is a ferromagnetic order when the lowest 
flat band is half-filled. Our results are consistent with 
this rigorous provement. In addition, in Ref . [l^ . quan- 
tum anomalous Hall effect in a flat-band ferromagnet on 
a two dimensional decorated lattice with spin-orbit cou- 
pling was also predicted. 

The paper is organized as follows. Firstly, we write 
down the Hamiltonian of the spinful Haldane model on 
honeycomb lattice with topological flat band in Sec. II. In 
Sec. HI, we do mean field calculations and get the phase 
diagram. In Sec. IV, we discuss magnetic properties and 
topological properties for the the half filling case. Then in 
Sec.V, we discuss the magnetic properties and the topo- 
logical properties for the the quarter filling case. Next 
we discuss the other doping cases in VI. Finally, the con- 
clusions are given in Sec. VII. 



II. THE MODEL HAMILTONIAN 
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In this paper we start from the spinful Haldane model 
with the on-site Coulomb interaction with the Hamilto- 
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FIG. 1: (Color online) The illustration of the honeycomb lat- 
tice. The red and black dots represent the bipartite lattice. 
The green lines, dash blue lines and dash magenta lines are 
the nearest hopping, t, the next nearest hopping t' and the 
next next nearest hopping, t" , respectively. 

nian as[l,[il,[3 

where U is the on-site Coulomb repulsive interaction 
strength and /i is the chemical potential, t, t' and t" are 
the nearest neighbor, the next nearest neighbor and the 
next next nearest neighbor hoppings, respectively (See 
Fig.l). We all know that the honeycomb lattice is a bi- 
partite lattice which is represent by the red and black 
dots in Fig.l, and the green lines, dash blue lines and 
dash magenta lines are the nearest hoppings, the next 
nearest hoppings and the next next nearest hoppings, 
respectively. In the flat-band limit we have t' = 0.6<, 
t" = —0.58^0. (pij = 0-4:TT is a complex phase to the next 
nearest neighbor hopping, of which the positive phase is 
set to be clockwiseQ. 

For the free electrons case (U/t = 0), the lower energy 
band becomes flat with nonzero Chcrn numbers Q ~ 2 
and the flatness ratio is about 50 (See Fig. 2). From Fig. 2, 
we can see that (a) is the totally energy bands, spin up 
and spin down are degenerate due to the spin-rotation 
symmetry and (b) is the amplification of the flat band. 
The flat band still has dispersion but the bandwidth is 
much smaller than the band gap. Fig. 3 illustrates the ex- 
istence of edge states in this Q — 2 topological insulator. 
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FIG. 2: (Color online) The free electronic dispersion {U/t = 
0) at t' = 0.6t, t" = -0.58t and (j>ij = OAn. (a) is the two 
bands, of which each band is doubly degenerate, (b) shows 
the details of the flat band. One can see that the bandwidth 
is much smaller than the band gap. 




FIG. 3: (Color online) The Zigzag edge states along x- 
direction at U/t = 0, t' = 0.6t, t" = -0.58t and (j>ij = OAtt. 
The black dash line is the position of chemical potential. 



III. MEAN FIELD APPROACH 

Next, we consider the interacting case oi U/t ^ 0. 
An interesting question is : what is the ground state for 
the spinful Haldane model with the on-site interaction in 
the TFB limit when the U/t ^ 0? As discussed above, 
people know that for the case of an interacting fermion 
model with flat-band, the ground state is always a FM 
order. On the other hand, near the half filling, the super- 
exchange effect will lead to an antiferromagnetic (AF) 
spin interaction between two local spin moments nearby. 
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Thus to answer this question, we will study the compe- 
tition between AF order and FM order by the following 
mean-field ansatzs : 



nian can be written as: 



1 



and 



(cUc,,,) = -[n + (-l)VM^;^] 



(2) 



(3) 



where Map is the staggered antiferromagnetic order pa- 
rameter and Mp is the ferromagnetic order parameter, n 
is the average particle number on each site. For the case 
of spin- up and spin-down, we have a = 1 and cr = — 1, 
respectively. 



A. AF order 

Firstly we consider the instability of AF order. For the 
AF order, we have Map 7^ 0, now the effective Hamilto- 

I 



(4) 



where Aa = UMap/^ and tr^ is Pauli matrix, n = 
1 — d is the average particle number on each site and d 
is the hole concentration. After Fourier transformation, 
the Hamiltonian in momentum space becomes 



H 



k 



where ^' 



alt ^it ^i-i '^U ' ^'^^ 



(5) 



hAk 



/-t'iC~D)+AA a- \ 

Ck -t'{C + D)~AA 

^t'{C-D)^AA 6 

V -t' {C + D)+ Aa J 



V 



where = -i + iB) - t" (E + iF). 

Then we can get the spectrum of the electrons as 



EA,k,i,±^-t'C+J\^k\+{t'D±A 



where 



EA,k,2,±=-t'C- ^J\Ck\ +{t'D±AAr 

A = 2 cos {kx/2) cos (^\/3ky/2^ + cos 
i? = 2 sin {kx/2) cos ( Viky/2 j — sin kx 



C = 2 cos (/) ^ cos (k • hi) 

3 

D ~ 2 sin (/) sin (k • b^) 

E = 2 cos fc^ cos + cos 2fcj^ 

_F = —2 sin fcj; cos (^y^ki^ + sin 2/03; 

where is the next nearest neighbor vectors. 



(6) 



(7) 



At different doping, by minimizing the free energy of 
the mean field, the self-consistent equations at zero tem- 
perature in the reduced BZ areflTl [isj: 



NsMap 



l-d 



E 

EA,k,2,~<fi, 

2A^., 



E 

t'D-A 



t'D + Aj 



2i 



kl 



E 



E 1] 



(8) 
(9) 



where ^^^.^ — 1^ — ^ and Ng is the number of the unit 
cells and 



&2 



\^k\' + {t'D + AA)\ 



+ {t'D~ AAf. 



In Fig. 4, we give the electronic dispersions for AF order 
at U/t = 7.0 for half-filling. Comparing with Fig. 2 and 
Fig. 4, we can see that AF order parameter can change 
the flat band into a dispersive one. 
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B. FM order 




Secondly we consider the instability of FM order. For 
FM order, we have Mp 7^ 0, now the effective Hamilto- 
nian can be written as 



J2 ^Fcja.c, + Y^(^-fA clc,, + h.c. (10) 



where Af = UAIf/2. After Fourier transformation, the 
Hamiltonian in momentum space becomes 



FIG. 4: (Color online) The electronic dispersions for AF order 
at U/t = 7.0, t' = 0.6t, t" = -0.58t and = OAtt for half 
filling. We can see that the AF order parameter changes the 
fiat band into a dispersive one. 



where 



H 



(11) 



J 



-t' (C-D)- Af a 

Ck -t'{C + D)-AF 







\ 



-t' {C-D) + Af a 

Ck -t'{C + D) + AF) 



Then we can get the spectrum of the electrons in FM 
order as 



EF,kx±^-t'C~AF±^M+{t'DY (12 



In Fig. 5, we give the electronic dispersions for FM or- 
der at J7/t = 7.0 for half- filling. Comparing with Fig. 2 
and Fig. 5, we can see that FM order parameter only splits 
the energy of the flat band without changing the flatness 
ratio. 



EF,k,2,± = -t'C + Af± V 161 + {t'Dy (13) 

At different doping, by minimizing the free energy of 
the mean field, the self-consistent equations at zero tem- 
perature in the reduced BZ are 



1 



1 



NsMf ^ 2 

+ E ^ E (14) 
+ E 1+ E 1] (15) 

where /x^^,^ = — ^^ and Ng is the number of the unit 
cells. 



C. Global phase diagram 

Then by using mean field approach above we can 
get the global phase diagram for different doping in 
Fig. 6 when the parameters are t' = Q.6t, t" = — 0.58i, 
~ 0.47r. From Fig. 6, we can see that there exist seven 
phases: normal metal, FM-metal, AF-metal, Q = 1 FM 
ordered topological Insulator, Q = 2 topological insu- 
lator with quantized anomalous Hall effect (we may also 
call it QAH state). A- type antiferromagnetic-spin density 
wave (A-type AF-SDW) topological insulator and trivial 
AF order. 

At half filling case d = 0, in the weakly coupling limit 
{U/t < 5.932), the ground state is a Q = 2 topological 
insulator with flat-band. In the strong coupling region 
{U/t > 5.932), due to Maf 7^ 0, the ground state be- 
comes an AF-SDW insulator. And at half filling case, 
the AF order always has lower energy than FM order. 
When after doping, the system is away from half filling. 
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FIG. 5: (Color online) The electronic dispersions for FM order 
at U/t = 7.0, t' = 0.6t, t" = -0.58t and 0ij = 0.47r for half 
filling. We can see that FM order parameter splits the two 
degenerate fiat bands without changing their flat ratio. 




FIG. 7: (Color online) The ground state energies of every 
site for AF order (Green line), FM order (Red line) and on 
spin order (black) per site at half filling for the parameters as 
t' = 0.6f, t" = -0.58t <f),j = 0.47r. 
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situation becomes much different : the FM order is really 
a Q = 1 topological insulator. For other hole concentra- 
tion cases d 7^ 0.5, d 7^ 1, the ground state becomes a 
FM metal. When the interaction strength becomes big- 
ger, FM order is stable in a large region in the global 
phase diagram due to the flatness of the fermion disper- 
sion. Only near half filling, the ground state energy for 
every site of AF order may be lower than FM order at 
some intermediate interaction strength (See below). And 
there exist first order quantum phase transitions between 
FM order and AF order (See the red dot line in Fig. 6). 

In the following, we will give a detailed discussion on 
the cases of d = 0, c? = 0.5 and 0, 0.5, respectively. 



IV. HALF FILLING CASE : d = 



FIG. 6: (Color online) The global phase diagram for the pa- 
rameters to be t' = 0.6t, t" = -0.58t and <f)ij = 0.47r. At half 
filling case d = 0, there are three phases: Q = 2 topological 
insulator with quantized anomalous Hall effect, A-TSDW and 
trivial AF order. For the quarter filling case; d = 0.5, there 
are two phases : paramagnetic Metal, Q = 1 topological Insu- 
lator with FM order. For other cases, there are three phases 
for the region of < d < 0.5: paramagnetic Metal, FM-metal 
and AF-metal while there are two phases for the region of 
0.5 < d < 1: paramagnetic Metal and FM-metal. 



d 7^ 0. Now the repulsive interaction is highly relevant 
for a fermion system with flat-band. And small inter- 
action will lead to an FM order. Thus in the phase 
diagram, away from half filling (d ^ 0), there exists a 
quantum phase transition between paramagnetic metal- 
lic phase and FM order for little interaction strength due 
to a big flatness ratio. And we flnd a narrow paramag- 
netic metallic phase. At hole concentration d = 0.5, the 



Firstly we consider the half filling case. At half filling 
d = 0, the particle number on average site is n = 1 and 
the chemical potential ^ = U /2 {l^^ff = O) . From Fig. 6, 
we can see that at d = the ground state is always 
dominated by the AF-SDW, for the ground state energy 
of AF order per site is lower than the FM order (See 
Fig. 7). In the weakly interaction limit {U/t < 5.932), 
the ground state is a Q = 2 topological insulator with 
flat-band. In the strong interaction region {U/t > 5.932), 
due to Map 0, the ground state becomes an AF-SDW 
insulator. 

Next we plot a phase diagram at half filling with dif- 
ferent values of t" /t in Fig. 8 (the purple dash line de- 
notes t" = —0.58t, the dispersion of largest flat ra- 
tio). In Fig. 8, there are three quantum phase tran- 
sitions {U/t)pM, {U/t)AF, and a topological quantum 
phase transition (the red line) that divide four phases: 
Q = 2 topological insulator, FM order, A-type topolog- 
ical AF-SDW (A-TSDW) and trivial AF order. For the 
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FIG. 8: (Color online) The phase diagram at half filling: there 
are four phases, QAH, FM-SDW, A-TSDW and Trivial AF 
order. The dash blue line represent the FM phase transition 
{U/t)FM, the black line represent the AF phase transition 
{U/t)AF, the red line represent a topological quantum phase 
transition. The purple dash line is the case with largest flat 
ratio, t' = 0.6f, t" = -0.58t, <?!>ij = QA-k. 



flat-band case (the purple dash line of t" ~ — 0.58t), one 
can see that the Q = 2 topological insulator occurs in the 
weak interaction, and with the increasing of U /i, firstly 
the system turns into A-TSDW state with Maf 0, 
and then turns into the AF insulator with trivial topo- 
logical properties. Now there only exist two types of 
phase transitions: one is the magnetic phase transition 
(the black line in Fig. 8) between a magnetic order state 
with Maf and a non- magnetic state with Maf = 0, 
the other is a topological quantum phase transition (the 
red line in Fig. 8) which is characterized by the condi- 
tion of zero fermion energy gaps (see the blue dot in 
Fig. 6). And in Fig. 9 and Fig. 10, we plot the staggered 
magnetization and the energy gap for the case with the 
parameters t' = 0.6t, t" — — 0.58t and (pij = QAtt at half 
filling. 

In Fig. 8, there is a dash blue line which is the boundary 
between Q — 2 topological insulator and FM order. At 
half filling case, the ground state energy of AF order 
per site is lower than the FM order. While away from 
t" = — 0.58i (the dispersion of largest flat ratio), there 
may still exist a narrow window of an FM order for t" /t > 
-0.55. 

At the half filling case, when the AF order parameter 
is non-zero, the electronic band structure will be reor- 
ganized. Thus in AF order, the spectrums of electrons 
are dispersive. In addition, we go beyond the mean field 
theory to study the collective spin fluctuations by ran- 
dom phase approximation (RPA) method. We obtain 
the dispersion of spin excitations for the ground state. 
See Fig.ll. From Fig.ll, one can see that in the AF or- 



FIG. 9: (Color online) The staggered magnetization Maf for 
the half filling case with the parameters t' = 0.6t, t" = — 0.58t 
and — OAtt. There are two quantum phase transitions (the 
black dash lines) : one is the magnetic phase transition at 
(U/t)AF, the other is a topological quantum phase transition 
with gap closing. 
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FIG. 10: (Color online) The energy gap of electrons for 
the half filling case with the parameters to be t' = 0.6t, 
t" = —0.58t and (pij — OAtt. There are two quantum phase 
transitions (the black dash lines) : one is the magnetic phase 
transition at {U/t)AF, the other is a topological quantum 
phase transition with gap closing. 



der where U/t > 5.932, {U/t = 9.0 in trivial AF-SDW 
order and U/t ~ 7.0 in A-TSDW order) there exist gap- 
Icss Goldstone modes - spin waves. We also calculate the 
dispersions of the spin wave of the AF order away from 
the flat-band case, t' = 0.6t, t" = -QMt, = OAn 
with U/t = 7.0 (blue line) and U/t = 9.0 (pink line). 
From this results, we can say that the spin waves are 
dispersive and have linear behavior, a;(fc) oc q around F 
pont, and also don't show anomalous behaviors at the 
flat band case {t" = -0.58t) in both trivial AF-SDW 
order and A-TSDW order. 
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FIG. 11: (Color online) The dispersions of the spin wave of 
the AF order for the case of largest flat ratio of the parameters 
are t' = 0.6t, t" = -0.58t, 4>ij = QAtv with U /t = 7.0 (balck 
line) and U /t = 9.0 (red line) at half filling. The dispersions 
of the spin wave of the AF order away from the flat-band 
case are also given : t' — 0.6t, t" = — 0.54f, (jtij = QAn with 
U/t = 7.0 (blue line) and U/t = 9.0 (pink line). 



FIG. 13; (Color online) The electronic enegy gap at the quar- 
ter filling (d — 0.5) for the case with the parameters t' — 0.6t, 
t" — — 0.58t and (j>ij = OAtt. There have two phases: the nor- 
mal paramagnetic metal (the Green region) and the Q = 1 
topological insulator with FM order. 




FIG. 12: (Color online) The staggered magnetization Mf at 
the quarter filling (d — 0.5) for the case with the parameters 
t' = O.Qt, t" = -0.58t and = OAtt. There have two 
phases: the normal paramagnetic metal (the Green region) 
and the Q — 1 topological insulator with FM order. In FM 
order, the ferromagnetization saturates suddenly. 



V. QUARTER FILLING CASE : d = 0.5 

Secondly we consider the quarter filling case: d = 0.5. 
In the global phase diagram in Fig. 6, there exists an im- 
portant line - the d = 0.5 line. We found that at the quar- 
ter filling the ground state is a Q = 1 topological insula- 
tor induced by the FM order Mp ^ 0. In Fig. 6 we can see 
that there exist two phases divided by a quantum phase 
transition aX U/t = 0.1 for d = 0.5 : the paramagnetic 



FIG. 14: (Color online) The electronic dispersion at the quar- 
ter filling (d — 0.5) for the case with the parameters t' = 0.6t, 
t" = -0.58t and ^tj = OAtv, U/t = 5.0. The chemical poten- 
tial fief f (the green plane) locates at the middle of the energy 
gap. The ground state is a Q = 1 topological insulator with 
filled lowest topological flat band. 

metallic state and the Q = 1 topological insulator. In the 
weakly interacting case {U/t < 0.1), the ground state is 
paramagnetic metallic state with flat-band. Small inter- 
action {U/t > 0.1) will drive the ground state to a Q = 1 
topological insulator with FM order. From this result 
one may guess that in flat-band limit (the flatness ratio 
is infinite), there exists only Q = 1 topological insulator 
for d = 0.5 and quantum phase transition between the 
paramagnetic metallic state and the Q = I topological 
insulator occurs at U/t = 0+ . This result is consistent to 
the rigorous proof in Ref. fl3| , in which it was proved that 
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FIG. 15: (Color online) The Zigzag edge state along x- 
direction at U /t = 5.0 with t' = QM, t" = -0.58t, = OAtt 
for the quarter filling case d = 0.5. The black dash line is the 
position of chemical potential. 
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FIG. 17: (Color online) The ground state energies of AF order 
and FM order per site for the case of d = 0.3 at t' = 0.6t, 
t" = -0.58t and = OAtt. The red line is the energy of FM 
order and the black line is the energy of AF order. The blue 
spots A and B denote the quantum phase transitions between 
AF order and FM order. 

firms the topological property of this Q = 1 topological 
insulator. In this phase, we also go beyond the mean 
field theory to study the collective spin fluctuations by 
RPA method (See Fig. 16). In Fig. 16, we can see that 
although the dispersion of electrons is flat, the collective 
spin modes are dispersive due the dispersive of the Mat- 
subara Green's function (See the detailed calculation in 
Appendix). However, at low energy limit, the situation 
of FM order is much different from that of the half filling 
case of AF-SDW order parameter : around the F point, 
the dispersion shows quadratic behavior instead of a lin- 
ear one as uj{k) oc q^. 



FIG. 16: (Color online) The dispersions of the spin wave of 
the FM order for the parameters are t' = 0.6t, t" — — 0.58f, 
(jjij = OAn with U/t = 3.0 and U/t = 6.0 at d = 0.5. 



the ground state of the Hubbard model with topological 
flat-bands is a ferromagnetic order when the lowest flat 
band is half-filled (that corresponds to the quarter filling 
case here). 

In Fig. 12 and Fig. 13, we show the ferromagnetization 
Mp and the energy gap for the case of c? = 0.5. In 
Q = 1 topological insulator, the ferromagnetization sud- 
denly jumps to 0.5 and gets saturated. In Fig. 14, we plot 
the electronic dispersion of the case d = 0.5 in Q = 1 
topological insulator region. From Fig. 14, we can see 
clearly that the chemical potential n^^^ located at the 
middle of the energy gap, and the lowest band is filled, 
while other bands are empty. We also calculate the edge 
states of Q = 1 topological insulator. Fig. 15 illustrates 
the edge states of the topological state. This result con- 



VI. OTHER CASES 

Thirdly, we consider other cases, d 7^ 0, d 7^ 0.5. Here 
we take d = 0.3 as an example. 

From Fig. 6, we can see that in the weak interaction 
case {U/t < 0.13), the ground state is normal metallic 
state with flat-band. Small interaction {U/t > 0.13) will 
drive the ground state to a FM metallic state. The phase 
transition between normal metal with Mp = and FM- 
metal with Mp 7^ is first order. The FM order param- 
eter jumps up at U/t = 0.13. When we increase inter- 
action strength, the ground state energy for every site of 
the FM order is larger than AF order, the ground state 
becomes the AF order crossing a first order phase transi- 
tion at U/t = 4.73. An interesting property is that when 
we further increase interaction strength the ground state 
energy for every site of the FM order is lower than AF or- 
der again, the ground state re-enter the FM order cross- 
ing another first order phase transition at U/t = 8.55. 
In Fig. 17, we give the ground state energy comparation 
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VII. CONCLUSION AND DISCUSSION 
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FIG. 18: ((Color online) The magnetization {Mf or Maf) for 
the doping d = 0.3 at t' = 0.6t, t" = -0.58i and 0ij = 0.47r. 
There have three phases: paramagnetic Metal (the Green re- 
gion), FM-metal and AF-metal. 




In this paper, we investigated the magnetic properties 
(including fcrromagnctism and antiferromagnctism) of a 
correlated topological insulator with flat band - the inter- 
acting spinful Haldanc model by the mean field approach. 
We found that the AF order is more stable than FM or- 
der at (or near) half-filling due to energy gaining from 
super-exchange effect; while away from half- filling, FM 
order is more stable than AF order due to the flat band. 
In particular, at doping d — 0.5, the system becomes a 
Q ~ 1 topological insulator which is induced by the FM 
order. In addition, we found that in FM order the degen- 
erate flat bands split without changing the flatness ratio, 
while in AF-SDW order, the electronic band structure 
changes totally and there doesn't exist flat-band at all. 
We also used the random-phase-approximation to calcu- 
late the dispersion of spin collective modes and found 
that in both FM order and AF order, the spin collective 
modes are dispersive. In AF order, the spin waves show 
linear dispersion while in FM order, the spin waves show 
quadratic dispersion. 
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(a) FM case: U/t=2.0 d=0.3 




(b) AF case: U/t=6.0 d=0.3 



FIG. 19: ((Color online) The electronic dispersion of the FM 
order and AF order at the doping d — 0.3 for t' — 0.6t, 
t" = -0.58t and (f)ij = 0.47r. (a) is the FM case and (b) is the 
AF case, we can see clearly that ^e// (the green plane) locate 
in the spectrum. 



between AF order and FM order per site at zero tempera- 
ture when d = 0.3, and we can see that there indeed have 
two transitions points A and B. In Fig. 18 we show the 
magnetic order parameters. In Fig. 19 we also plot the 
electronic dispersion at the same doping to compare the 
different dispersion in AF order and that in FM order. 
Correspondingly, the electronic band structure changes 
when phase transitions occur. 



VIII. APPENDIX: RPA CALCULATIONS FOR 
SPIN COLLECTIVE MODES 



In this appendix, we use random-phase-approximation 
(RPA) to calculate the spin dispersion at doping d ~ Q 
and d = 0.5. The spin excitations are obtained from 
the poles of the transverse spin susceptibility tensor, x, 
which is defined, in Matsubara form, as 

xtj {q,i^n)=l dre'-'^'-^xti (9''^)' 
Jo 

where 

and i,j = a,b label the two sublattices and = 



E 



4 



^fe ^fc + ,,t'^f=,4. ' ^-q Efc '-fc„, I '-fc,t r- 

erators for each sublattice. Using the Wick's theorem, 
we can derive the leading order spin suspectibility in the 
AF phase and FM phase as 



\ ^ denote the spin op- 



+-(0) / • \ 



xG7(fc- 



q.iujyi + iujrn) ■ (16) 



10 



In AF phase, the Matsubara Green's functions are : 

1 / 1 I t'_D+o-A 

2 1 



tibiUty tensor in AF order from the Dyson equation 



icu — EcrJ {k) 



■ 1 — "u",? 



3=± 



2C.,i(fe) 



iuj — Ea-j (k) ' 



Ck 



(17) 

(18) 
(19) 
(20) 



(21) 



where / denotes the 2x2 identity matrix. The poles of 
the spin susceptibihty tensor, corresponding to the spin 
excitations, are then obtained from the condition : 



Det / - C/x° 



0. 



(22) 



where 



E^^± (fc) = ±J {tA + t"Ef + [tB + t"Ff + {t'D + al\Af 



t'C - /ie// 



We note that the tensorial nature of the spin susceptibil- 
ity is a consequence of two sublattices here. 

The dispersion of spin collective modes with frequency 
(jj and momentum q is determined by above condition 
(j22p after performing the analytic continuation iu) — >■ 
Lo + zO~^. Finally the spin wave dispersion determinant 
is given by 



e.,± (fc) ^ ± V {tA + t"Ef + (tB + t"Ff + {t'D + al^A? 



= -{tA + t"E) -i{tB + t"F) 
We may use the RPA approach to get the spin suscep- 



i^+_(<7,c.) = l-C/x+-(°)-C/x, 



-(0) 



+ ^ Xta ^ '^Xbb ^ ^ Xab ^ 'Xfca 

= 0, 



-(0) 



(23) 



where the elements of the spin susceptibility tensor 
X+-(°) are : 



_ t'D{k)+AA \ f, I t'D{k-q)-AA \ f-, _ t'D{k)+AA \ _ t'D{k-q)-AA 
-(0)/, ..N 1 V ^tW Uik-q) ) V i,{k) A g^Cfe-g) J \n( p (,.^^ 

N.^iy 4(^„-i?,,+(fc-g)+i?t,-(fc)) + 4(..„-i?;,_(fc-g)+i?,,_(fc)) ""^'-^^^^ 

i(fc)+AA\ _ t'D{k-q)-AA \ (. _ t'-D(fc)+A.4 ^ (. _ t'D(k^q)-l 



k 



^ ^ t'D{k)+AA \^ t'D(k-q)-AA \^ t'D{k)+AA ^ t'D(k^q)-AA ' 

4 {ojn - E^^- {k-q) + St,+ (fc)) ^ 4 (t^„ - Ei^- {k-q) + ^t,- (fc)) 



'(-i?^,_(fc-g)) 

(24) 



+-(0) ^ \ 
Xfcb {q,^n) = 




, t'D{k)+AA \ _ t'D{k~q)~AA \ {, I t'n(fc)+A.4 ^ , t'D{k~q)~AA \ 



4 (cj„ - (fc - g) + i;t,- (fc)) 



4(cj„-S;,_(fc-g) + i;t,_(fc)) 



0(-£;t,-(fc)) 



' /-I _ t'D{k)+AA \ , t^D(fc-g)-AA \ /-, f'D(fc)+AA ^ , t' D(k~q) 
4(L^„-S^,_(fc-q) + St,+ (fc)) 



4(t^„-£;^,_(fc-(7) + St,-(fc)) 




11 



k — q 



TV, ^4et(fc)Ci(fc-g) 



1 



1 



1 



LJn - -B;,- (fc - g) + ^t,- (fc) - {k-q)+ E^^^ (k) ^ 

1 



w„ - (fc - q) + E^^_ (fc) a;„ - [k - q) + E^^+ (fc) 

I 



ei-Ei,^ik-q)) 



(-i^t,- 
(26) 



TV, ^ 4^t (fc) (fc - 9) L - Ei.- {k-q) + St,- (fc) cj„ - E^^+ (k - q) + E^,^ [k) , 



1 



1 



w„ - (fc - g) + £;t,- (fc) a;„ - E^,^ {k - q) + St,+ (fc) 

I 



6* (-i^t,- 
(27) 



In FM order, we may use similar approach to calculate tive modes, 
the spin suspcctibility and the dispersion of spin coUcc- 
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